Non-equilibrium Fluctuation Relations in a Quantum Coherent Conductor 
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We experimentaUy demonstrate the vahdity of non-equihbrium fluctuation relations by using 
a quantum coherent conductor. In equilibrium the fluctuation-dissipation relation leads to the 
correlation between current and current noise at the conductor, namely, Johnson-Nyquist relation. 
When the conductor is voltage-biased so that the non-linear regime is entered, the fluctuation 
theorem has predicted similar non-equilibrium fluctuation relations, which hold true even when the 
Onsager-Casmir relations are broken in magnetic fields. Our experiments qualitatively validate the 
predictions as the first evidence of this theorem in the non-equilibrium quantum regime. 

PACS numbers: 05.40.-a, 72.70.+m, 73.23.-b, 85.35.Ds 



The fluctuation-dissipation relation m is a central con- 
cept in physics. It was triggered by the description of 
Brownian motion by Einstein Q, who claimed that the 
response of a physical system to an external force is pro- 
portional to its equilibrium fluctuation. In electrical cir- 
cuits this fact manifests itself as the Johnson-Nyquist 
(JN) relation 0] in such a way that the conduction 
through a conductor is proportional to its current noise 
in equilibrium. The linear response theory Q founded 
on this relation together with the Onsager-Casimir recip- 
rocal relations provides a powerful tool to describe a 
variety of physical systems. This is, however, justified 
only when the systems are close to equilibrium, and the 
generalization to the non-equilibrium regime has been 
long sought. 

Generally the current (/) passing through a conductor 
can be expressed as a polynomial of the bias voltage iV) 
as follows; 



(1) 



where the first term represents Ohm's law with conduc- 
tance Gi (Such an expansion in V is assumed to be valid 
for mesoscopic transport Q). While in a mesoscopic con- 
ductor Gi is directly related to the transmission Q, the 
higher-order coefficients (G2, G3,...) convey information 
on electron-electron interactions in a voltage-biased con- 
ductor as shown recently [iMl- Similarly to Eq. (1), 



a polynomial of V for the current noise power S gener- 
ated in the conductor, namely the variance of I, can be 
expressed as; 



S ^ Sq + SiV + ^^32^ + 



(2) 



The coefficients of the first terms in Eqs. (1) and (2) are 
linked by the JN relation 6*0 = iksTGi, where ks is 
the Boltzmann constant and T is the temperature of the 
conductor 0, Now, an essential question arises: 
Are there also correlations between the higher orders (in 
voltage) in current / and noise power 5? 

Here we demonstrate a non-equilibrium fluctuation re- 
lation in a quantum coherent regime. We find a propor- 
tionality between Si and G2, which corresponds to the 
next order correlation beyond the JN relation between Sq 
and Gi. The present relation is asymmetric in magnetic 
fields unlike the Onsager-Gasimir symmetry and is valid 
even in the presence of interactions. Such a direct link 
between non-linearity in voltage and fluctuations out-of- 
equilibrium has been predicted theoretically 15-1^ by 
applying the fluctuation theorem [2^ to a quantum co- 
herent conductor. Our experiments qualitatively validate 
the predictions but disagree on the quantitative level. 

We used an Aharonov-Bohm (AB) ring as a typical co- 
herent conductor. Figure 1(a) shows an atomic force mi- 
croscope (AFM) image of the AB ring fabricated by local 
oxidation using an AFM [iH on a GaAs/AlGaAs het- 
erostructure two-dimensional electron gas (2DEG) (the 
electron density 3.7 x 10"^^ cm~^, the mobility 2.7 x 10^ 
cm^/Vs, and the electron mean free path 2.7 /ini at zero 
back-gate voltage) as well as the experimental setup for 
the two-terminal measurement in a dilution refrigerator. 
As our 2DEG has a back gate to tune the electron den- 
sity, the conductance of the AB ring can be modulated 
by the back gate voltage Vg and the magnetic fleld B 
(the AB effect). Figure 1(b) shows the conductance as 
a function of Vg and B, displaying clear AB oscillations 
with an oscillation period being 25 mT, in agreement 
with the ring radius of 230 nm f]3\. The conductance 
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FIG. 1. (a) AFM image of the AB ring with the DC and 
noise measurement setup in the dilution refrigerator whose 
base temperature is 45 mK. The in-plane gates defined by the 
oxide lines are grounded in this experiment, (b) Conductance 
of the AB ring as a function of Vg and B. (c) Gi (left axis) 
and So (right axis) as a function of Vg at B = T. (d) So 
is plotted as a function of Gi. The solid line indicates the 
JN relation of So = AksTGi with T = 125 mK. (e) G2 (left 
axis) and S\ (right axis) as a function of Vg. (f) Si is plotted 
as a function of G2. The solid line is the result of the fitting 
(Si = 3.64 X 4kBTG2 with T = 125 mK). 



of the ring at S = T and = V is 3.1 in units of 
e^/Zi (25.8 kn)-'^. The visibility of the AB oscillations 
around B = T is 0.13 (see Fig. 3(a)). 

In addition to the DC measurement, we performed a 
noise measurement as follows (see Fig. 1(a)). The volt- 
age fluctuation across the sample on the resonant circuit, 
whose resonance frequency is about 3.0 MHz with the 
bandwidth of ^ 140 kHz, is extracted as an output sig- 
nal of the cryogenic amplifier 22h24|| . The time-domain 
signal is then captured by a two-channel digitizer, and 
is converted to spectral density data via FFT, where the 
cross-correlation technique is adopted to increase the res- 
olution. By numerically fitting the obtained resonance 
peak, the current noise power spectral density S is de- 
duced as performed in Ref. '23']. 

We deduce the coefficients in Eqs. (1) and (2) by nu- 
merically fitting the measured current I and current noise 
power spectral density S as polynomials of V [2^. In 
the analysis we set the bias window to \eV\ < 50 /xeV, 
where Joule heating is negligible as confirmed in previous 
noise measurements [iJ] (Two examples of the analysis 
are shown in Fig. 2) . Figure 1 (c) shows the conductance 
Gi (left axis) and the equilibrium noise power density Sq 
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(right axis) as a function of at i? = T. Due to elec- 
tron interferences, Gi varies as Vg changes. The behavior 



FIG. 2. (a) (b) (Lower panels) Two examples of I-V charac- 
teristics obtained at different Vg's, where ^G2V^ components 
deduced from the numerical fitting are superposed as dashed 
curves. (Upper panels) The corresponding S, where the re- 
sults of the polynomial fitting and SiV components are su- 
perposed as solid and dashed curves, respectively. The non- 
equilibrium noises are smaller than the equilibrium ones at 
certain bias voltages indicated by the arrows. The signs of 
G2 and Si are correlated (negative in (a) and positive in (b)). 



of 5*0 perfectly follows that of Gi as expected from the JN 
relation. The proportionality between Gi and Sq shown 
in Fig. 1(d) indicates that So = iksTGi is fulfilled with 
an electron temperature of T = 125 mK. 

Now let us focus on the next order in voltage V. Fig- 
ure 1(e) shows G2 (left axis) and Si (right axis) as a 
function of Vg. Both G2 and Si vary reflecting elec- 
tron interferences. Moreover, their variations are cor- 
related (Fig. 1(f)) with a correlation factor being 0.88 
between the two. When the result is expressed in the 
form 5*1 = 4:kBTG2a" ("0" denotes the zero-field) based 
on the assumption that the functional form is similar to 
the JN relation, we obtain a" = 3.64 ± 0.54 by a nu- 
merical analysis 0, where the error bar indicates the 
95% confidence interval (see the solid line in Fig. 1(f)). 
While the magnitude of a° will be discussed later, two 
remarks are made. First, a proportionality between G2 
and Si is already contained in the scattering theory for 
non-interacting systems [1], which yields general results 
to express the current and its noise as a function of the 
bias voltage, the transmission of the conductor, and the 
Fermi distribution of the reservoirs at a given tempera- 
ture. However, we will show later that the present rela- 
tion is valid beyond the description based on this picture 
in terms of the breaking of the Onsager-Casimir symme- 
try. Second, the finite values of Si indicate that there is 
a certain bias voltage range where we measure "negative 
excess noise" , namely the non-equilibrium noise smaller 
than the equilibrium one [1, H, 27 1. Indeed, depending 
on the coefficients {Si, S'2,...) in Eq. (2), the value of the 
bias voltage to give the minimum noise power is shifted 
from zero by a value of the order of 10 /iV [T§| as shown 
in Fig. 2. 
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The Onsager-Casimir reciprocity is a fundamental re- 
lation to characterize the response of a system close to 
equilibrium Figure 3(a) shows Gi (left axis) and 
(right axis) as a function of i? at = 0.02 V. Both of 
them, displaying a periodic oscillations due to the AB 
effect, are symmetric with regard to i? = T, namely 
Gi{B) = Gi{-B) and SoiB) = So{-B). Furthermore, 
they are proportional to each other (Fig. 3(b)). This 
is the direct consequence of the JN relation and the 
Onsager-Casimir reciprocity. On the other hand, the 
next order coefficient G2 has no such symmetry as shown 
in Figs. 3(c) and 3(e) (blue curves), where the symmet- 
ric components Gf (B) = G2{B) -\- G2{—B) and anti- 
symmetric components Gf{B) = G2{B) — G2{—B) are 
presented, respectively. The emergence of a finite at 
B 7^ is in agreement with previous works on the break- 
ing of the Onsager-Casimir reciprocity in the non-linear 

[IMl- 



transport regime 



It was attributed to electron- 
electron interactions in mesoscopic conductors [lol . 
since, in the Landauer picture for non-interacting elec- 
trons [§] , the transmission probability of electrons always 
obeys the Onsager-Casimir reciprocity regardless of ap- 
plied bias voltages, leading to a vanishing anti-symmetric 
term G^. 

Here we show that the non-equilibrium relation be- 
tween G2 and Si is valid even when the Onsager-Casimir 
relation is broken. In Figs. 3(c) and 3(e) the symmet- 
ric components 5'f (i?) = Si{B) + Si{~B) and anti- 
symmetric components Si{B) = Si{B) — Si{—B) as 
a function of B are superposed in red curves, respec- 
tively. In addition to the proportionality between Gf 
and Si (see Fig. 3(d)), there exists a clear proportion- 
ality between Gf and Sf as shown in Fig. 3(f). The 
correlation factors are 0.94 and 0.88 for the symmetric 
and anti-symmetric parts, respectively. For the symmet- 
ric part, = 2.69j^Q'35 is obtained in the expression of 



Sf = AksTG^a^, while = 
the anti-symmetric relation S^ 



2.18I^f| is deduced in 
= AkBTGfa^. Impor- 
tantly, the anti-symmetric part is non-trivial in that it 
is the consequence of the departure of the system from 
the Onsager-Casimir symmetry, and a signature of non- 
equilibrium and non-linearity. Nevertheless, the higher- 
order correlation still exists there as well as in the sym- 
metric part, but with different coefficients. 

What is the reason for the correlations in the non- 
linear transport regime? Here we give an intuitive pic- 
ture of higher order correlations following recent theo- 
ries 15l4l9l|. Electron transport can be viewed as the 
electron exchange process between two reservoirs via a 
conductor. Consider the probability P{Q) that one reser- 
voir gains Q electrons from the other as a consequence 
of the exchange (in the absence of a magnetic field, for 
simplicity). As time reversal symmetry, particle-number 
conservation, and energy conservation are required over 
the whole system including the reservoirs, P{Q) should 
satisfy P{Q) = P(-g) exp(fi^Q) This means 
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FIG. 3. (a) Gi (left axis) and So (right axis) as a function 
of B at Vg ^ 0.02 V. (b) So as a function of Gi. The solid 
line So = AksTGi is the JN relation at T = 125 mK. (c) 
Symmetric components in the second-order coefficients Gf 
(left axis) and Sf (right axis) as a function of B. (d) Sf as 
a function of Gf . The solid line is the result of the fitting 
corresponding to Sf = 2.69 x iksTCi for T = 125 mK. (e) 
Anti-symmetric components in the second-order coefficients 
G2 (left axis) and Si (right axis) as a function of B. (f) Si 
as a function of G2*. The solid line is the result of the fitting 
corresponding to Si = 2.18 x AksTGi for T = 125 mK. 



that the probability of a Q-electrons exchange process 
and that of its time-reversed process (— Q-electrons ex- 
change) are linked with a factor difference of exp( j^Q), 
reflecting a probability difference of the initial states 
for the two exchange processes. Note that this equal- 
ity based on the micro-reversibility is valid even in the 
non-equilibrium situation and provides the basis of the 
fluctuation theorem [23|. As P{Q) contains all the infor- 
mation associated with electron transport through the 
conductor, the strong constraint posed by this equality 
on the cumulants of Q yields a number of correlations 
starting with the JN relation. For the second order, the 



equality predicts Si = 2kBTG2 [15l - ll9j . In the pres- 
ence of magnetic fields, the fiuctuation relations that do 
not rely on the microreversibility out of equilibrium were 
recently pointed out [l3, [3] , where the anti-symmetric 
relation is expressed by S^ - 2kBTGf = C^/SksT. 
Here, Cq is the anti-symmetric part of the third cumu- 
lant ("skewness") at equilibrium. With the restriction 
of micro-reversibility, this anti-symmetric relation [l6l | 
is reduced to Si ~ e/cs^G^ with another relation 
sf = Cq /2kBT. In both cases, the symmetric relation 
is given by = 2kBTG^. 
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FIG. 4. (a) Values of So/d obtained at 125, 300, and 
450 mK. The solid line shows the JN relation, (b) Values of 
Sf/Gi at 125, 300, and 450 mK. The solid line is the theoret- 
ical expectation (S'f /Gf = 2fcsr) [ll-[il. The dashed line is 
the result of the linear fitting that gives 

(c) Values of /G2 at 125, 300, and 450 mK. The solid line 
is the theoretical expectation (S^ /G2 = QknT) [l^.lT9l|. 



At zero-field, the observed linearity between Si and G2 
qualitatively agrees with the above discussion, although 
the theory does not well reproduce the obtained ratio 
Si/2kBTG2 = 2a° ~ 7. At finite magnetic field as well, 
while the expected proportionality is obtained, the ra- 
tio {Sf/2kBTGl = 2a^ - 5) is larger than in theory. 
The present observation for a° and is perfectly repro- 
ducible in our experiments performed in several differ- 
ent conditions for different magnetic fields and gate volt- 
ages. Concerning the temperature dependence, although 
the ratio Sq/G^ follows the JN relation (Fig. 4(a)), 
the ratio S'f/Gf shown in Fig. 4(b) depends linearly 
on temperature as expected, but the slope is similarly 
larger than the theoretical prediction. This implies the 
necessity of further theoretical work that takes realis- 
tic situations into account [2^. In contrast, the result 
for the anti-symmetric component {S^ /dhsTG^ ~ 1.5) 



is in better agreement with theory Si — GfesrCj in 
spite of a possible independent contribution of Cq when 
only the universal relation without micro- reversibility 
(S^ 



2kBTG^ 



C^/SfcsT) is assumed [I7|,ll8|. The [20 



agreement is also clear in the temperature dependence of 
the ratio S'j^/G^ as shown in Fig. 4(c). This observa- 
tion suggests that micro-reversibility in magnetic field is 
likely to be validated in the present experiment. 

In conclusion we have experimentally proven the pres- 
ence of the direct link between the non-linear response 
and the non-equilibrium fluctuation in the AB ring, as 
theoretically predicted on the basis of the fluctuation 
theorem. While our demonstration was performed for 
the simplest case in a normal coherent conductor with- 
out energy relaxation in the scattering region, the present 
finding will further open an applicability of the fluctua- 
tion theorem in the quantum coherent regime and in the 
presence of magnetic fields. 
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Appendix: Simple deduction of the higher order 
correlations 



We explain the fluctuation theorem using the simplest 
setup. We consider a mesoscopic conductor, say a quan- 
tum point contact, where the two quantum wires are cou- 
pled by tunneling. For simplicity, no magnetic field is ap- 
plied. The present system is described by the following 
Hamiltonian 



H = Hl + Hr + H, 



LR, 



(A.3) 



where and Hn are the Hamiltonian of the left and 
right quantum wires and Hlh is the tunneling part be- 
tween them. The initial density matrix is decoupled into 
the equilibrium states of each wire, where the left and 
right wires are assumed to have equal temperature 
and have chemical potentials hl and hr, respectively. 
Then the whole density matrix is described by 



Pinitial — / ^ PriL^nu 

\nL,nR){nL,nii\, (A. 4) 



where and Zr are the normalization factors, and 
|'^L,"-fl) defines the state that ul and ur electrons are 
present inside the left and right wires with the eigenen- 
ergies and i?„^ of Hl and Hr, respectively. The 
probability to find the state |n^, u'r) after a certain time 
T starting from the initial state \nL,nR) is expressed as 



(ni,nH)^(n'^,n'^) 



We note the time reversal symmetry 

We also note electron number conservation tll — n'j^ = 
— {ur — n^) and energy conservation satisfied at very 



large r: E^'^ - E^^ ~ -(^n'^ - Em,)- Using these time 
reversal symmetry and conservation laws, we find the 
relation 

P, > / / , — P, , : \ , ^pMriL-n'^) 

where A is an affinity A = (3{fiL — Pr)- The probabil- 
ity that transmitted electron number is Q, is defined as 
P{Q) = E„^,„„,„i,,„;,/'(«x,,n«)^(„'^,«'^)'5(Q-(ni-n^)). 
A direct consequence from the relation (|A.6P is "Fluctu- 
ation Theorem" : 

P{Q) = P{-Q)e^Q. (A.6) 
Now let us discuss the higher order correlations be- 
tween the current and its noise power, which are ad- 
dressed in the present paper. With Fluctuation Theorem 
(|A.6[) . we find the following identity 

= J dQQP{Q) = - J dQQP{Q)e-^^ 



-{Q) + MQ')-^{Q') + ' 



(A.7) 



Furthermore, we note that (Q") is a function of A, i.e., 

(Q") = (Q")o + + ^(Q")2 + • • • . (A.8) 

Comparing order by order with respect to A, we find 
infinite number of relationships among these quantities, 
some of which are given as 



(g')i = (g)2. 



(A.9) 
(A.IO) 



Average current / and current noise power S are defined 
as / {Q)/t and 5 = 2{{Q^) - {Q)'^)/t. The first re- 
lation (jA.9[) is equivalent to the fluctuation dissipation 

50 = AkBTGi [G. GaUavotti, Phys. Rev. Lett. 77, 
4334 - 4337 (1996)], and the second relation (|AlO|) is to 

51 = 2kBTG2 (See main text for the deflnitions of Gi, 
G2, Sq, and 5*1). More systematic and exact derivation 
for these relations including the finite magnetic field case 
is performed by using a cumulant generating function [K. 
Saito & Y. Utsumi, Phys. Rev. B 78, 115429 (2008)]. 



